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SUMMARY 


Using the cold neutron scattering technique about 300 phonons have been determined in a single 
aluminium crystal at room temperature to define 10 pairs of dispersion curves. Investigations 
have been made of the variation of frequencies, phonon line widths and multi-phonon spectra in 
the temperature range 293 < 7’ < 932°K. For a particular direction in the crystal lattice it is 
shown that the frequencies vary about 15 % over this temperature range. The line widths are of 
such a magnitude that the derived phonon mean free paths vary from about 5 phonon wave lengths 
at 600°K to about 1.5 phonon wave lengths at 930°K. The observed multi-phonon spectra are 
found to agree with calculated differential cross sections in the incoherent approximation. 


I. Introduction 


As indicated by theoretical work by Waller and Froman [1], and Placzek and 
Van Hove [2] slow neutrons should be a useful tool to investigate the properties of 
the phonons in solids. Later successful experiments by Brockhouse et al., and Hughes 
and Palevsky et al. have shown that the careful analysis of the energy and angular 
distribution of the slow neutrons scattered inelastically in fact makes it possible to 
study single phonon processes and determine the dispersion relations w = (q) in 
the case of coherent inelastic scattering [3] and to derive the vibrational frequency 
spectrum in the case of incoherent inelastic scattering [4]. In general the inelastic 
scattering processes are in the coherent case governed by the energy and momentum 
conservation laws 


<A €= +1 for energy gain 
Hout — Lin = Eh 2 or eee, ‘i fice (la) 
Kout — Kin = + 2 Gi (1 b) 
i=1 


where E =h? k?/2m represents the neutron energy, the q;:s are the phonon wave vectors 
belonging to the frequencies «;. As long as the temperature of the single crystal under 
study is lower than the Debye temperature, Op, of the crystal, the most pronounced 
inelastic coherent scattering effect is by far the single phonon scattering with m = 1. 
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At higher crystal temperatures, however, the multi-phonon scattering becomes in- 
creasingly important, the cross section for these processes increasing rapidly with 
temperature [5]. Disregarding a minor effect of singularities in the two-phonon scatter- 
ing [6] (m = 2) the neutron spectra expected from higher phonon scattering (m > 2) 
are continuous without any pronounced maxima and minima as in the case of one- 
phonon scattering. This should at least be true for the incoherent scattering. In the 
case of coherent inelastic scattering there might possibly occur some weak intensity 
maxima and minima in the scattering picture for the two-phonon case but for three- 
phonon and higher processes any sharp structure should be washed out by the averag- 
ing effect taking place in the scattering. 

The conclusion of these theoretical predictions is that there is a good chance that 
the technique of the inelastic coherent scattering of slow neutrons could be used even 
at high crystal temperatures, 7’ > 6p, to give information on the single phonons just 
because of the pronounced discontinuous nature of the one-phonon processes and 
because of the continuous nature of the multi-phonon scattering processes. 

With these ideas in mind a study of the lattice vibrations in an aluminium single 
crystal from room temperature, 293°K, up to the melting pint, 932°K, was under- 
taken, using the cold neutron scattering technique. The reason why aluminium was 
chosen as sample was that its cross section is all coherent, which is of particular 
importance if one wants to avoid complications in the interpretation of the continuous 
multi-phonon spectra obtained at higher sample temperatures. 


II. Experimental method 


a. General 


In general the experimental problem encountered in inelastic neutron scattering 
studies is twofold: first to produce a monochromatic incoming neutron beam to 
define k;, of equation (1) and second to perform an angular and energy analysis of 
the scattered neutrons defining the vectors k,,,. In the present experiment the cold 
neutron technique using a beryllium filter to produce a spectrum with the well known 
cut off at 3.95 A was used to define the ingoing neutron wave length [3, 4, 7, 8]. The 
disadvantage with this technique is the large width of this spectrum—the full width 
at half maximum about | A—which in general leaves the question open as to what 
the effective ingoing neutron energy actually is. Particularly this might be serious 
for small energy gain processes, where this gain is only of the order of (1-2)- Ein. 
The intensity transmitted by a filter is, however, very high, about 1% of the total 
thermal neutron flux in a reactor beam being transmitted if the filter is cooled to 
liquid nitrogen temperature [7]. Using a special technique of interpretation of data 
it 1s possible to overcome the difficulty presented by the width of the spectrum as 
will be discussed in the next section. 

Using the filter technique it is also possible to work with large samples, a fact 
which of course improves the scattered intensity. In the present experiment a beam 
area of 5 x 10 cm? was utilized and the single aluminium crystal used had a diameter 
of 6.5 cm and a length of ~12 em. This is an advantage as far as study of single 
phonon processes is concerned. However, in the study of the broad multi-phonon 
scattered spectra, the multiple scattering effect in big samples might be of such a 
magnitude that the use of smaller samples is unavoidable. 

The complete experimental set up involving the filter, chopper and time-of-flight 
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Fig. 1. Spectrometer set up at the Stockholm reactor. 
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analyzer used at the Stockholm reactor has been described in detail elsewhere [9, 10], 
so it should here only be emphasized that the arrangement permits the scattered spec- 
tra to be observed at various angles of observation between 0° and 90° with respect 
to the incoming beam. A photograph showing some essential parts of the spectro- 
meter is given in figure 1. 


b. Interpretation of scattered spectra 


As just mentioned the main problem using the filter technique is to associate the 
right ingoing wave length with any particular scattered intensity maximum. To 
simplify the generally rather complicated shape of the scattered spectra all the ob- 
servations to be described in the following were made in a principal plane, the 001 
plane, of the aluminium single crystal, i.e. the crystal was oriented with its [001] 
direction horizontally, coinciding with the axis of rotation of the spectrometer arm. 
This arrangement simplifies the observations because one of the “‘transverse’’ 
branches (more correctly: one of the low frequency branches) is associated with zero 
scattering intensity, leaving us in general with two intensity maxima for any direc- 
tion of observation and any crystal orientation as long as the observation is made in 
this principal plane. The observed branches then are one high and one low frequency 
branch corresponding to vibrations of longitudinal and transverse nature respectively. 
Only in directions of high symmetry these vibrations are pure. A typical scattering 
picture obtained at room temperature for a 90° angle of observation is shown in 
figure 2. As shown in this figure the ingoing broad spectrum might be reproduced in 
two ways by one-phonon scattering (the multi-phonon scattering is at room tempera- 
ture associated with a very low cross section) either (a) all the scattered intensity is 
concentrated in one narrow peak, (b) the scattered spectrum very nearly reproduces 
the shape of the ingoing cold neutron spectrum including the sharp cut off. In certain 
directions other intermediate and complicated alternatives might occur, but these 
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Fig. 3. Explanation of the spectrum given in figure 2. a. Narrow maximum formation. b. Broad 
maximum formation. Ingoing spectrum given along the vector k,,. 


are the principal shapes of the scattered spectra. To explain these two important 
scattering pictures the scattering surfaces generated for different k,, vectors are 
given in figure 3.1 In general each incoming neutron energy or wave vector k{? 
corresponds to a particular scattered neutron energy or wave vector k{},. If the 
direction of observation, gy, is kept fix and if observations are made for different 
crystal orientations 6, each of the k&}, vectors will generate a scattering surface S. 
If the scattering surfaces so generated happen to be nearly parallel to the direction 
of ingoing wave vectors k{? it is seen that all the scattered wave vectors k{), very 
nearly have the same length (case a). We have a focusing effect, which permits us 
to use for instance the longest k,, vector corresponding to Aj, = 3.952 A to define 
the energy gain as indicated in figure 2 for the intensity maximum (a). Also it is 
obvious that this mechanism of formation of an intensity maximum is equivalent to 
a considerable intensity gain as all the available intensity of the ingoing spectrum 
is used in the formation of one single peak. On the other hand if the scattering surfaces 
generated by the different k{} vectors make a larger angle with the vectors kf? the 
vectors k‘), vary in lengths as indicated in figure 3b. If now the differential cross 


u 
section for the scattering process is such that 


4/7 tt 


6 (kin —>kout) > 0 (Kin —>Kout) > 0 (Kin >kout) and so on, 


1 In figure 3 actually parts of the intersections of the different scattering surfaces with the 
principal plane under discussion are given. 
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Fig. 4. Experimentally determined scattering surfaces at two different angles of observation and 
room temperature. 


the scattered spectrum will reproduce the form of the ingoing beryllium filtered 
spectrum as shown. This permits us to use the wave length shift of the beryluum break 
at 3.952 A to define the energy gain as indicated in figure 2 for the intensity maximum 
(b). This discussion may of course be carried on for other possible scattering pictures, 
and it is found that with the exception of directions near singular directions [11] it is 
always possible to use the wave length 3.952 A to define k,, and associate it with the 
corresponding scattered wave length defining k,,,. To carry such an investigation 
through parallel to the collection of experimental data, it is necessary to have some 
idea of the frequency spectrum of the crystal or rather to know the isofrequency 
surfaces to a first approximation. Fortunately these had been given by theoretical 
calculations [12] for the case of aluminium partly based on X-ray measurements [13], 
so these data were used in combination with the known cross section formulas [1.11] 
to analyze the observed intensity distributions. As described elsewhere (reference 14, 
the appendix) small corrections amounting to 1-3 % of the observed shifted beryllium 


cut off positions on the time-of-flight scale had sometimes to be made as an effect of 
the finite resolution of the instrument. 


c. Collection of data 


Knowing the principle of interpretation of the data the next experimental problem 
using a fix ingoing wave length, 3.952 A, is to fill the reciprocal lattice plane, 001, with 
experimentally determined phonon wave vectors q and corresponding frequencies w 
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to be able to construct the dispersion relations @ = co(q). This is, however, possible 
using the full ability of the experimental equipment. If the angle of observation 
increases, the momentum transfer in the scattering process, k,,,—k, =x, also 
increases for the vector x having a fix direction in the reciprocal lattice. This property 
of the scattering process immediately suggests the method to use: 


1. for a fix angle of observation, g, the crystal orientation is varied in small angular 
steps—in the actual observation in steps of 2.5°—and the scattering surfaces are 
determined for this angle of observation. 

. the same series of observations is repeated for different angles of observations, 
the scattering surfaces being determined for each angle q. 


bo 


In practice because of the symmetry properties of the crystal the observations can 
be limited to a 45° interval in the first Brillouin zone. Most of the observations to be 
described were made in the zone round the 220 reciprocal lattice point. The reason 
for this choice was that both the high and low frequency branches here give scattered 
neutron intensities concentrated into rather narrow peaks because of the focusing 
effect illustrated in figure 3a. Typical results gained at room temperature are illus- 
trated in figure 4. The observations at ~ = 90° were made over a wide region of 
angular settings of the crystal to test the ideas outlined above. The observations at 
g — 40° are typical for the series of measurements made at other angles m and similar 
observations were made at g = 83, 80, 70, 60, 50 and 30°. 


III. Room temperature results (T = 293° K) 


In the way described in the preceding section about 400 phonons were determined 
at room temperature of the crystal, all these phonons falling in the 001 plane. Be- 
cause of the symmetry of the crystal all of them are moved into a 45° sector of the 
first Brillouin zone falling between the [100] and [110] directions. With a few excep- 
tions for cases of weak intensity or for cases where the interpretation of the data was 
more complex, the accuracy of the measurements is such that the q vectors are 
determined to +(3-5)% and similarly the frequencies w also are known to an ac- 
curacy of + (3-5) % in general. With all the experimentally determined points moved 
into the 45° sector, lines were drawn from the 220 reciprocal lattice point with 5° 
intervals, in all 10 directions between [100] and [110] and including those, i.e. [100], 
[100] + 5°, [100] + 10° ... [110]. All experimental points falling within +1° from 
these lines were used to determine points on the 10 dispersion curves. In those cases 
when the experimentally determined phonon vectors did not fall within +1° from 
the selected direction the intersection between the scattering surface and this direc- 
tion was used to define the phonon on the dispersion curve. This resulted in about 
300 phonons being determined and the resulting m (q) relations are given in the table of 
appendix I. Examples of the results for three directions are given in figure 5. For 
the [100] and [110] directions the results of X-ray measurements [13] are also given 
for comparison. There seems to be a definite discrepancy for the high frequency, 
the ‘longitudinal’, branch, the neutron data giving shgtly higher frequencies. For 
the low energy branch, the “‘transverse”’, there is good agreement in the [100] direc- 
tion and a marked disagreement in the [110] direction, again the neturon data giving 
higher frequencies. ‘ 

Although no serious attempt was made to compare the present data to theory— 
mainly because it was felt that the fitting of trial models has little value until a new 
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Fig. 5. Dispersion curves determined in the present experiment compared to X-ray and theoretical 
data. 


fundamental treatment of the dynamical problem is made—a simple model was 
used to calculate the dispersion curves in the directions for which experimental 
dispersion relations had been found. The model used in this calculation! was that due 
to de Launey [15]. In this model central forces are supposed to act between an atom 
and its nearest neighbours corresponding to one force constant and between the 
atom and its next nearest neighbours corresponding to a second force constant. 
Furthermore the effect of the electron gas surrounding the ions is taken into account 
in such a way that the “impure” high, medium and low frequency vibration vectors 
are divided up into three new “‘pure’’ vectors of which one is longitudinal and two 
are transverse. The electron gas is supposed to resist compression, i.e. to affect the 
pure longitudinal components of the waves but not the transverse components. A 
correction to the forces is thus introduced using the electron gas bulk modulus. The 
dynamical problem is completely determined by the three elastic constants ¢,,, Cys 
and ¢,,. Even if the physical meaning of such a model did not seem quite clear, for 
our purposes such a model appeared to have certain advantages as the variation 
of the frequencies with temperature could be determined from the known temperature 
dependence of the elastic constants [16] in a first approximation. The results of these 
calculations are shown as solid curves in figure 5. The values of the elastic constants 
at 293°K were used. Although there is no good quantitative agreement between theory 


and experiment the tendency of this relatively simple theory is in amazingly good 
agreement with experiment. 


* Carried out by Dr. Alf Sjélander, now at Argonne National Laboratory, U.S.A. 
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IV. Results gained at elevated temperatures 293 < T < 932° K 
a. Experimental technique 


In making neutron scattering experiments at elevated temperatures one is con- 
fronted with the problem of constructing a reliable heating system capable to main- 
tain the crystal at a high constant temperature without gradients and yet this oven 
should not scatter or absorb neutrons to any appreciable extent. In a first version, 
for use in the temperature range 300 < 7’ < 800°K, a thin mica sheet, ~0.5 mm thick, 
was wrapped round the crystal as insulator and then a flat resistance band was 
wired directly on the mica. The crystal with this simple direct heating system was 
placed in a polished, larger steel container of 0.5 mm wall thickness closed on all 
sides to secure a certain degree of heat insulation from the surroundings. The 
temperature was controlled by three thermocouples penetrating to different depths 
into the crystal. The crystal was mounted between two pins in such a way that it 
could be rotated round its [001] direction. 

In a later version when the crystal was heated up to the melting point a thin steel 
container was made fitting exactly over the crystal to hold and contain it and again 
the same technique was used to insulate and heat it. The thermocouples used to 
measure and control the temperature by automatic regulation of the heating current 
ensured a temperature definition to +5°K all over the crystal volume. 

Even if great care was taken to make the oven thin relative to neutron scattering, 
some background scattering was observed particularly for the highest temperatures. 
Therefore separate background runs with an empty oven system were made for each 
temperature used and for different angles of observation as the most important 
background caused by the oven was the inelastic scattering from iron in it. This 
background never exceeded 15-20% of the measured effect from the aluminium 
crystal. 


b. Results 


As it could be expected the multi-phonon—particularly two-phonon—scattering 
events should get an increasing importance at higher crystal temperatures, it was 
necessary to make a careful choice of the position in the reciprocal lattice where one 
should try to study the one-phonon process. In accordance with the discussion given 
in section IIb, it would be of great advantage if the intensity of the ingoing spectrum 
could be concentrated into one narrow intensity maximum in the one-phonon scatter- 
ing process as illustrated in figure 3a. Fortunately this happens for the longitudinal 
as well as for the transverse branches in a region of the reciprocal lattice near the 
220 reciprocal lattice point. Therefore a first series of runs was made at an angle of 
observation of g = 90° and a crystal orientation defined by the angle 6 ~ 65° and first 
for crystal temperatures in the range 293 < 7 <793°K. The results of these first 
measurements are given in figure 6. In this figure are also given results of a later 
run at 7 =910°K and a final run taken at the melting point a 932°K. From this 
first series of runs two effects were clearly observed and the possibility of a third 


effect appeared: 


1. A multi-phonon spectrum rapidly grows up for increasing temperature and be- 


comes dominant at about 700°K. 
2. The one-phonon lines shift their position in such a way that the frequency , 


defined by the neutron energy gain, decreases with temperature. 
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3, Also the one-phonon lines seemed to broaden. At least this effect was observed 
at the highest temperatures. This effect was, however, difficult to estimate from 
this first run because of the strong multi-phonon scattering. 


From these first observations it was thus obvious that the change of the dispersion 
curves ® =«(q) with temperature could be measured, if the direction of observa- 
tions was chosen in the proper way. Also the multiphonon spectrum as well as the 
line width could be studied, if one could find a position in the reciprocal lattice 
where the scattered neutron intensity caused by the “longitudinal” phonon was 
zero. Hf furthermore the frequency of the phonon studied could be chosen relatively 
low this would mean that the multi-phonon spectrum and the one-phonon line 
would be separated from each other on the wave length or time-of-flight scale. 
Actually it was possible to find such a place in the reciprocal lattice choosing a direc- 
tion of observation of g = 50° and a crystal orientation defined by 6 ~ 60°. As the 
phonon line widths and line shifts are intimately connected to the question of an- 
harmonicity of the forces acting between atoms, we went into a detailed study of 
the temperature effects. The observation of the three phenomena will now be treated 
separately. 


1. Variation of frequencies w with temperature 


Knowing the dispersion relations at room temperature and with the information 
from the first runs at higher temperatures (figure 6) it was possible to choose such 
combinations of angles of observations, g, and crystal orientations, 6, that for each 
temperature the phonons could be determined along one direction in the reciprocal 
lattice, i.e. the dispersion curve could be determined directly for each temperature 
investigated. The direction selected was one making an angle of ~18° with the [100] 
direction in the 001 plane corresponding to the [3810] direction. The phonon wave 
vectors determined in this way always fell within +2° from the desired direction. 
Thus the dispersion relations were determined for this direction for temperatures 
293, 387, 486, 700 and 793°K for the high frequency branch and for 293 and 793°K 
for the low frequency branch. Figure 7 shows these dispersion curves as well as the 
ones calculated using the model of the Launey. The temperature variation of the 
theoretical curves was derived from the measured temperature dependence of the 
elastic constants [16]. In calculating the experimental curves from the measured 
scattered neutron spectra due consideration was given to the volume expansion of the 
crystal (contraction of the reciprocal lattice). In determining the position on the 
time-of-flight scale of the one-phonon peaks, it was carefully studied, if this position 
could be affected by the multi-phonon spectrum on which the peaks were super- 
imposed (compare figure 6). The true shape of the multi-phonon spectrum could be 
nea to good accuracy using the theory of Sjélander (compare following sec- 
tions). 

Again it is striking how the simple model of de Launey predicts a variation of 
frequency which gives the general tendency reasonably well. Quantitatively there is 
no agreement as the predicted temperature variation for a fixed q/dmax is about 50% 
higher than the measured. 

To investigate how the frequency varies with temperature for a constant q/max; 
the frequency values read from the experimental dispersion curves drawn in figure 7 
for different constant values of ¢/Gmax were plotted against temperature. The result 
is given in figure 8 where also are added the result of some studies at higher tempera- 
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Fig. 7. Experimentally determined dispersion curves in the direction [310] at higher temperatures 
(left) compared to the theoretical ones (right). 
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tures for ¢/ max = 0.3. The frequencies vary linearly with temperature even up to the 
highest temperatures near the melting point. Also the slope of the lines is independent 
of the value of ¢/dmax in the region, 0.3 <q/qmax <0.8, we have investigated. The 
frequency variation with temperature is well described by the formula 


w (q, 7’) = @ (q, To) — (9.5+ 1.0): 109 (7 — T,) (2) 


So far this temperature variation of the frequencies is thus studied for one direction 
in the 001 plane of the reciprocal lattice. To attempt to estimate the agreement or 
disagreement of these results with the measured temperature variation of the Debye 
temperature, 9p, of aluminium, let us assume that the frequencies for any direction 
in the crystal vary in the same way. If we take a frequency of about 4-10! rad/sec 
as an average frequency in the spectrum, we find the percentage variation of this 
frequency to be about 12% in going from 300°K to 800°K. According to X-ray 
measurements of Owen and Williams [17] the Debye temperature varies about 15 % 
over the same temperature range. The present neutron results thus do not seem to 
contradict X-ray results. 


2. Variation of line width with temperature 

To be able to study the line width variation it was necessary to find a low frequency 
phonon, which in the scattering process should produce an intensity maximum as 
narrow as possible and as well separated from the multi-phonon spectrum as possible. 
As mentioned earlier in this section it was possible to find such a combination of 
angles g and @ that these conditions were reasonably well fulfilled. Figure 9 shows 
the scattering pictures obtained for m = 50° and 6 = 60° and for the temperatures 293, 
486, 700, 854 and 923°K. All backgrounds are subtracted. Similar observations were 
made at a large number of temperatures between 773 and 932°K. A striking effect is 
that the definition of the one-phonon lines persists even up to the immediate neigh- 
bourhood of the melting point. The multi-phonon spectrum is reasonably well sepa- 
rated from the single one-phonon maximum, such that the multi-phonon part may be 
subtracted from the intensity under the one-phonon peak as indicated in the figure 
by the dotted lines at the base of the one-phonon peaks. To perform this subtraction 
the theoretical multi-phonon spectra discussed in the next section were of great assis- 
tance. After this subtraction the widhts of the one-phonon line were analyzed in the 
following way: 

The line width observed is built up of three different contributions: a. the width of 
the line caused by the wave length spread in the ingoing spectrum 0/, (compare 
figure 3b),! b. the resolution width, 6/,, which was 90 ys or 0.116 A, c. the natural line 
width 6/,, which we want to determine. As there are three factors which make up 
the total, observed line width 6A), it is a reasonable approximation to use the formula 


5 An = VO 15 — 61 — 67; (3) 
to calculate 6/,, i.e. to assume the components to have approximately Gaussian 


shape. For temperatures below 600°K the line width 6/, obviously was so small, that 


it could not be observed, i.e. 6A) = V ba? a5 O72. Of the two contributions to dA,, the 
resolution width, 6/,, certainly stays constant. It remains, however, to show that 


1 Actually the ‘‘one’”’ phonon peak is not just formed by one single phonon but by a band of 


q , 
phonons corresponding to a 15-20 % spread, from —! ~0.55 to ~ 0.65. 
Imax Imax 
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Fig. 9. Types of observed scattered spectra used to analyze line widths. 
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Fig. 10. Phonon line width as a function of temperature. 


0A; also stays constant, such that any observed increase in width really is caused by 
an increase in the natural line width, 6/,. A variation of 6/; should mean that the 
angle between the vectors k$}, and the scattering surfaces S should change with 
temperature (see figure 3a or b). By differentation of the energy and momentum 
relations (1) and using the measured dispersion curves it is, however, possible to show 
that the change in the line width caused by this effect certainly is less than 4°K over 
the temperature range of our measurements. As already mentioned, for temperatures 
below 600°K the width of the line stayed constant and equal to 135 ys or 0.174 A. 
The accuracy to which this could be determined was about +5 ws. From 600°K and 
on a monotonous increase could be observed ending up with a width of about 180 jus 
or 0.232 A near the melting point. The analyzed results with the phonon line width 
expressed in °K are given in figure 10 and appendix IT. Included are also a few results, 
which could be gained by analyzing in the same way the width of a “longitudinal” 
phonon of similar frequency. 

In order to estimate the mean life time and the mean free path of the phonons 
studied, a calculation of the line shape assuming an exponential damping was carried 
out.! Assuming the phonon waves described by 


i(qa-R-w.t —y,|t 
U,—ce ™ ORO Ys ltl (4) 


and working the differential scattering cross section out using this assumption, it 

turns out that the cross section instead of 6-functions of the form 6(@ — a, — € as) 

for the case y, = 0 will now contain a factor describing the line shape 
1 

OO Ses)" : 


Vs 


(5) 


iy els 


1 
S(O —O,— €@;)= : 
a ( 


where fi (@ — wy) is the neutron energy gain (€ = + 1), im, is the energy of the phonon. 
This result may be interpreted in the way that when damping is present the energy 


1 The calculation carried out by Dr. Alf Sjélander. 
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Fig. 11. Phonon mean life time and mean free path as a function of temperature. 


conservation law is not strictly fulfilled. The line gets a width corresponding to energy 
exchange between phonons during the damping. It is thus obvious that the full width, 
AE, at half maximum for this line shape is 2y;:h. The mean life time, Af, is from 
eq. (4) I/ys such that AH At = 2h. Using the measured line widths, 6 7°K, it is then 
easy to calculate the mean life time At =2h/k;6 7’, where kz is the Boltzmann constant. 
Knowing the mean life time, this is multiplied by the wave velocity to get mean free 
path. Because the phonons studied are both to be found on the almost linear portion 
of the dispersion curves—q/qmax ~ 0.6 for the ‘‘transverse” phonon with @ ~ 2.4- 1018 
rad/sec and q/qmax ~ 0.3 for the ‘longitudinal’? phonon witha ~ 3.0- 1013 rad/sec—it is 
immaterial if the phase or group velocity is used to calculate the mean free path. 
The mean life times and mean free paths thus derived are given in figure 11. As the 
wave lengths of the transverse and longitudinal phonons studied are about 6 A and 
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144A respectively, these observations indicate that the mean free paths vary from 
about 5 times the phonon wave length at 600°K to about 1.5 times the phonon wave 
length near the melting point. The damping of the phonon waves is thus very strong, 
indicating the harmonic approximation to be rather doubtful at these crystal tem- 
peratures. 


3. The multi-phonon spectrum 


As already illustrated best by figure 9 the multi-phonon spectrum rapidly grows 
up when the crystal temperature exceeds the Debye temperature of the crystal. 
A calculation of this spectrum for the coherent scattering case is of almost prohibitive 
complexity. The cross section has, however, been given in convenient form in the 
incoherent approximation. Sjélander [5] has given the differential cross section in 
the form of a phonon expansion (see ref. 5, formulas I1.28 and I1.29), which is just 
a form suitable for comparison with the present experiment. 

However, before the data may be compared to theory all backgrounds have to 
be carefully subtracted. In measuring a continuous, broad spectrum of this type using 
a crystal of large dimensions the problem of multiple scattering is of greatest im- 
portance. Therefore this phenomenon was studied at length both experimentally 
and theoretically. As indicated by calculations the multiply scattered neutron spectra 
should be of a shape similar to the multi-phonon spectrum giving proportionally 
greater contributions to the observed spectrum for smaller neutron energy gains. To 
investigate the multiple effect experimentally several different experiments were 
tried. In one case an aluminium poly-crystal of the same size as the single crystal was 
cut into two halves along its cylinder axis. At temperatures of 600°K and 800°K 
first an upper half was irradiated by the now well collimated beam, this beam only 
seeing this upper half and the scattered spectrum was recorded at angles of observa- 
tion of g = 90° and 50°. Then the lower half was added but not irradiated directly 
and the scattered spectra were again recorded. When the spectra thus recorded 
were compared for constant gy, they were found to coincide within the statistical 
accuracy of about 3-4 %, so there was no observable change in intensity or spectrum 
shape caused by multiply scattered neutrons from the lower half. Then the same sort 
of experiment was repeated now reversed, i.e. with the lower half irradiated and the 
two different kinds of spectra were recorded. When compared they were again found 
to have the same shape within +4 % but the intensity for the case when both halves 
were present was found to have decreased by about 10%. This decrease is, however, 
what can be expected from the thickness and cross section, if the upper half in 
this case merely acts as an absorber and outscatterer for neutrons coming from the 
lower half. Thus there was no observable effect due to multiple scattering. 

In another experiment we irradiated again at 600°K and 800°K firstly a solid 
aluminium poly-crystal of cylinder shape and secondly an aluminium tube of the 
same diameter, 6 cm, as the solid cylinder. As the tube wall thickness was 0.6 cm 
in the case of the hollow cylinder, the multiple scattering effect should be much smaller 
than in the case of the solid cylinder. Again, however, the two scattered spectra as 
compared for constant angle of observation were of the same shape within the statis- 
tical accuracy. ; 

From all these negative results we could only draw the conclusions that either 
the multiple scattering is small or that the multiply scattered spectrum should have 
the same shape as the single scattered spectrum. Consequently we could not apply 
any corrections for the multiple scattering to the observed multi-phonon spectrum. 
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Fig. 12. Observed multi-phonon differential cross sections compared to the theoretical ones. 
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final uncertainty, however, is about 25%. The cross sections were also calculated 
absolutely for different crystal temperatures using Sjélander’s formulas. As the 
Debye temperature, 9p, enters critically into these calculations a linear variation 
of 0p from 6p = 380°K at T = 293°K to 6p = 325°K at T' = 923°K was assumed. This 
amounted to a 14% variation of p over this temperature range, a reasonable varia- 
tion in the light of these measurements and of the measurements of Owen and 
Williams. For comparison a calculation of the cross section at 7’ = 923°K was also 
made assuming 6p constant = 380°K. In calculating the cross section at the highest 
temperatures up to five-phonon processes were included. Also the cross section was 
integrated over the ingoing spectrum shape to be comparable with experimental 
data. Finally the experimental resolution function was applied to the theoretical 
curves. The result of these calculations are given in figure 12 as well as the experi- 
mental points. The experimental values are normalized to the theoretical at 7’ = 923°K 
as there was some disagreement (although it happened to be small) between the 
absolutely calculated and measured values. The cross section is given on an absolute 
energy (w) scale such that the ingoing spectrum is to be found below wm ~ 0.8: 10% 
rad/sec. As seen from this figure there is a rather satisfactory agreement between 
theory and measurement. The fact that the incoherent approximation so well de- 
scribes the cross section indicates that the coherence effects in two- and certainly in 
higher-phonon processes are averaged out in these multi-phonon scattering events. 
Also the curve shown for the case of 923°K and 6p = 380°K shows that the variation 
of 6p with temperature certainly has to be taken into account to explain the variation 
of the inelastic scattering cross section with temperature, a fact already found for 
total cross section measurements [18]. 
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APPENDIX I 


Frequencies for given phonon wave vectors. 


Direction in the ie 
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(a = 0° =[100}- 
direction, | oe 1028 Ox lOsss 
a = 45° =[110]- Gnas Padicee Imax rad/sec 
direction) 
0° 0.120 + 0.024 1.13 +0.05 0.330 + 0.017 1.53 £0.04 
Ge = oo LAS 0.160 + 0.030 1.55 = 0.07 0.333 + 0.017 1.85 + 0.05 
0.200 + 0.031 2.20 £0.10 0.432 + 0.017 2.13+0.09 
0.361 + 0.026 3.49 +0.12 0.598 + 0.017 2.87+0.10 
0.400 + 0.030 Ay zis O22 O6U7 sa 0017 2.89 70.11 
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0.881 + - 0.027 
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direction, OoalOm 
a= 45° =(110]- 4 Gmax rad/sec 
direction) 
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0.465 + 0.015 
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1.48 +0.05 
1.97 + 0.04 
2.44+0.11 
2.55 + 0.07 
2.52 +0.06 
2.48 + 0.06 
2.64 + 0.06 
2.84+0.07 
2.97+0.07 
2.97+0.18 
2.92+0.12 
2.94 +0.12 
3.45 +0.11 
3.53 + 0.16 
3.53 £0.11 


0.97 £0.03 
1.17 + 0.04 
1.234 0.04 
1.05 + 0.03 
1.40 + 0.04 
2.54 70.07 
2.49 + 0.06 
2.32 + 0.08 
2.82 + 0.08 
2.79 +£0.14 
2.78 + 0.08 
2.98 20.12 
3.26 + 0.13 
3.27 0.08 
Se One OLS 


1.02 + 0.03 
1.04 + 0.03 
1.24+0.04 
1.32 + 0.04 
1.39 + 0.04 
2.51 + 0.10 
2.49 40.12 
2.88 + 0.13 
2.83 + 0.14 
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ARKIV FOR FYSIK. Bd 17 nr 21 


Direction in the 
001-plane. Longitudinal Transverse 
(a =0° =[100]- 
. eee @ x 10-13 @ x 10-18 
aoe es rad/sec 4 Fmax rad/sec 
40° 0.773 + 0.039 5.29 + 0.22 0.570 + 0.016 2.64 + 0.06 
(ee NGo2vAS 0.778 + 0.031 5.25 + 0.19 0.594 + 0.016 2.57+0.14 
0.803 + 0.036 5.61 + 0.29 0.730 + 0.016 3.09 + 0.12 
0.864 + 0.027 5.42 + 0.26 0.734 + 0.016 3.07 + 0.07 
0.950 + 0.025 5.31 + 0.20 0.761 + 0.016 3.07 + 0.07 
0.961 = 0.028 5.16 + 0.23 0.931 + 0.016 3.43 + 0.07 
0.972 + 0.044 5.12 + 0.25 
45° 0.165 + 0.031 1.86 = 0.09 0.213 + 0.016 1.02 + 0.03 
Umax = 1-640 A-? 0.310 + 0.024 3.82 + 0.12 0.227 + 0.016 1.28 + 0.04 
0.445 + 0.031 4.67 0.18 0.483 + 0.019 2.48 + 0.10 
0.555 + 0.028 4.95 + 0.19 0.536 + 0.016 2.57 + 0.06 
0.643 + 0.028 5.24 40.18 0.539 = 0.022 2.83 + 0.09 
0.784 + 0.031 5.22 + 0.18 0.687 + 0.016 2.92 + 0.08 
0.843 + 0.028 5.50 + 0.25 0.845 + 0.016 3.23 £0.07 
1.000 = 0.028 5.31 40.21 
1.000 = 0.038 5.10 £0.23 
1.000 = 0.044 5.08 + 0.23 


APPENDIX II 


Observed line widths 6f,,, corresponding to 6A, in the text, and derived phonon 
line widths, 67, given in °K as a function of crystal temperature. 
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Longitudinal vibration 
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